
STATISTICAL TEST METHOD FOR OBJECTIVE VERIFICATION OF 
AUDITORY STEADY-STATE RESPONSES (ASSR) IN THE 

FREQUENCY DOMAIN 


2^ The invention relates to the field of objective measuring of hearing 
ability by using evoked auditory steady-state responses (ASSR). The 
proposed methods for objective verification of ASSR in the frequency 
domain can be employed in an ASSR-based hearing screening of a 
newborn as well as in objective audible threshold measuring based on 
ASSR. 

There are two different types of ASSR known: 

1. Click-evoked ASSR 

2. ASSR evoked through a continuous amplitude or frequency- 
modulated tone, which are also referred to as amplitude-modulation 
following response (AMFR). 

Both types of ASSR are described (Figure 1) in the frequency 
domain by so-called harmonics (one fundamental wave and several 
harmonic waves). The frequency of the fundamental wave corresponds to 
the click stimulus rate or the frequency of the modulation signal. The 
frequencies of the harmonic waves are multiples of the frequency of the 
fundamental wave. This means that the entire response is represented 
by a few spectral lines. The substantial portion of the noise power caused 
by the asynchronous electroencephalogram (asynchronous EEG) is 
concentrated, in contrast, on the spectral lines lying between the 
harmonics. 
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The objective verification of the ASSR occurs almost exclusively in 
the frequency domain. Each spectral line in the frequency domain is 
defined by a spectral amplitude and a phase angle. Several statistical 
methods are known for response verification in the spectral domain, 
which interpret in so-called one-sample tests only the phase angle or also 
the phase together with the amplitude of an individual spectral line, 
preferably the first harmonic (fundamental wave) (Stapells DR, Makeig S, 
Galambos R. Auditory steady-state responses: Threshold prediction using 
phase coherence. Electroencephalography and Clinical Neurophysiology 
1987;67:260-270; Valdes JL, Perez-Abalo MC, Martin V, Savio G, Sierra 
C, Rodriguez E, Lins O. Comparison of statistical indicators for the 
automatic detection of 80 Hz auditory steady state response (AMFR). Ear 
and Hearing 1997;18:420-429). The recorded time signal is transformed 
by epochs into the frequency domain for this purpose. The length of the 
transformed epochs must be selected in such a manner that the epoch is 
an exact integral multiple of the length of the period of the click stimulus 
rate or of the modulation frequency. The spectral line (fundamental 
frequency) corresponding to the click stimulus rate or the modulation 
frequency is searched and tested in the frequency spectrum that exists 
after the transformation. The advantage of the response verification in 
the frequency domain, compared to the direct verification in the time 
domain, is that the portion of the noise power in the spectral domain 
(represented by the spectral lines lying between the harmonics) does not 
interfere with the response detection since these spectral lines are not 
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included in the testing. The disadvantage of the so-called "one-sample 
tests" employed in the above-mentioned publication exists in the 
limitation of the statistical testing on the fundamental frequency. 
However, ASSR are not represented, in the rule, alone by the click- 
repetition frequency or the fundamental frequency corresponding to the 
modulation frequency, but they are also represented by one or several 
harmonic waves for which there is allotted a portion of the response 
signal power, which is not be ignored. An objective verification method 
limited only to the fundamental wave is therefore not optimal. 

Statistical testing methods working in the frequency domain are 
also described by Sturzebecher and Sturzebecher et al. (Sturzebecher E: 
"Method for hearing screening of newborn by means of steady-state 
response evoked with high click rate", European patent application EP 
01610060.4.; Sturzebecher, E, Cebulla M, Baag M, Thie R: "Verfahren zur 
objectiven frequenzspezifischen Hdrschwellenbestimmung mittels der 
Amplitude-Modulation Following Response (AMFR)", European patent 
application EP1099408 A2). So-called q-sample test are employed as 
statistical tests, which include the fundamental wave as well as the 
relevant harmonic waves for the statistical verification of the ASSR. We 
are dealing thereby with the "q-sample uniform scores test" known in the 
literature (Mardia KV. Statistics and directional data. Academic Press 
London and New York 1972) and a modification of this test proposed by 
Sturzebecher et al. (Sturzebecher, E, Cebulla M, Baag M, Thie R: 
"Verfahren zur objectiven frequenzspezifischen Hdrschwellenbestimmung 


mittels der Amplitude-Modulation Following Response (AMFR)", European 
patent application EP1099408 A2). 

The calculation rule for the "q-sample uniform scores test" 
described my Mardia, 1972, is: 

Let {Xj k ; 1< i < m, 1 < k < q} be a collection of random variables (phase 
angles cp ik ) ; q is the number of samples (spectral lines) with the sample 
size m (number of epochs), i.e. there are q x m = n phase angle values. 
The n phase values were ranked in a single sequence. Let r ik , i = 1, 
m, be the ranks of the phase angles in the /cth sample. 
The phase angles <Pii<are then replaced by the uniform scores 

n 

The test statistics used is 

m m 

with Ck =^J3 ik - St-Y&P* 
/=i i=i 

where 

r ik are the ranks of the n phase angles (n = q x m), q is the number of 
samples (number of included spectral lines) 

and m is the sample size (number of epochs). 

W is distributed as Chi-square with 2 (q - 1) degrees of freedom. 

As it can be seen from the calculation rule, only phase angles are 
used and the spectral amplitudes are not considered. Another loss of 


information is added: only ranks of phase angles are used in the 
calculation of the test value and not the phase angles themselves. This 
has, nevertheless, the advantage that the method is nonparametric; 
however, the result of the information loss taken in exchange is a lower 
test power. However, the test power should be as high as possible for 
hearing screening and for the objective audible threshold determination. 
Sturzebecher et al. have therefore devised a modification of the test 
(referred here as Test Modification 1), which considers also the spectral 
amplitudes in the form of ranks of amplitudes in addition to the phases: 
Test Modification 1 

Additionally to the phase angles, the spectral amplitudes A ik were 
taken into account. Like the phase angles, the spectral amplitudes Aj k are 
ranked in a single sequence: Let a ik , i = 1,..., m be the ranks of the 
spectral amplitude A ik in the kth sample. The phase angles cp ik were 

2-n: -r ik 

replaced by the uniform scores A* = — - — 

The test statistics used for the modified q-sample uniform scores test is 



with C* k = ]T a& • cos fc. . «S* = 2 



a ik -smfi ik 


and 


A ik = 


n 


where 
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r jk are the ranks of the n phase angles (n = q x m), q is the number of 
samples 

m is the sample size of the q samples and 

aik are the ranks of the n spectral amplitudes A jk . 

This already known test modification gains more information 
through inclusion of the spectral amplitudes than the "q-sample uniform 
scores test" of Mardia (1972) and it is still parameter-free because of the 
use of ranks instead of the real phase and amplitude values; however, 
working with ranks of phases and amplitudes instead of the actual values 
still means that the existing information is not used completely, which has 
the consequence that the test power is not optimal. t, 

It is the object of the invention to create additional diverse 
modifications of the known u q-sample uniform scores test" in which the 
available information of phases or the phases together with the spectral 
amplitudes are fully used. 

The following novel modifications of the known "q-sample uniform 
scores test" are proposed: 

Test Modification 2 

Only the phase angles are used; however, in contrast to the known 
"q-sample uniform scores test" of Mardia, 1972, one does not work with 
the ranks but with the phase angles computed through Fourier 
transformation. 


Let {Xix; 1< i < m, 1 < k < q} be a collection of random variables (phase 
angles <p jk ) ; q is the number of samples (spectral lines) with the sample 
size m (number of epochs), i.e. there are q x m = n phase angle values. 
The test statistics used is 


where 

q is the number of samples (number of includes spectral lines) 
and m is the sample size (number of epochs). 
Test Modification 3 

Spectral amplitudes and phase angles are used in the known 
modification; however, the phase angles are not ranked while the ranks 
for the spectral amplitudes are still entered in the test. 
The spectral amplitudes Aik are ranked in a single sequence; 
Let aik, i = 1,-, m be the ranks of the spectral amplitudes A ik in the kth 
sample. 

The test statistics used is 


W 


m 


m 


with Q=]Tcos^; 


s k =2] sin Pa 


7=1 



with 



m 


where 


q is the number of samples 
m is the sample size of the q samples and 
aik are the ranks of the n spectral amplitudes A ik 
<p ik are the phase angles. 
Test Modification 4 

In this case, there are used directly (unranked) the values of the 
phase angles computed by means of Fourier transformation and the 
spectral amplitudes as well. 
The test statistics used is 


where 

q is the number of samples 

m is the sample size of the q samples and 

A jk are the spectral amplitudes 

<Pik are the phase angles. 

A substantial problem in the use of the inventive test modification is 
the fact that the respective associated density function of the test values 
(probability density function) is unknown for the null hypothesis as a 
result of the performed modification. The critical test values necessary 
for the use of the test can therefore not be taken from the tables in the 



m 

with c;=/£4, 


ik ■ cos <Pik l = J A* ■ sin <p a 


m 


7=1 


current literature (the test values are necessary for the decision of 
positive or negative test results). 

The Monte Carlo simulation offers a known possibility for 
computation of the density function of the null hypothesis. A very large 
amount of pairs of random numbers are created hereby with a random 
number generator. A spectral amplitude and a phase angle are calculated 
from each pair of numbers and the statistical test is used thereafter. The 
distribution of the null hypothesis is calculated from the resulting large 
amount of test values. The searched critical test value can be read from 
the distribution. 

However, the normal distribution of spectral amplitudes and phases 
accepted in the simulation cannot be assumed in the case of real spectral 
amplitudes and phase angles. 

An element of the invention is therefore the following method for 
computation of the distribution of null hypothesis which takes into account 
the real distribution of spectral amplitudes and phase angles: 

Assumed is the presence of a great number (>100) of ASSR 
recordings whereby raw data of the derived electroencephalogram (EEG) 
has been continuously stored on the hard disk (approximately 200 epochs 
whereby the length of one epoch is approximately 1 second). As it is 
depicted in Figure 1, the response is limited to a few spectral lines (the 
fundamental wave is 160 Hz and the higher harmonic is a multiple of 160 
Hz). Since the spectral resolution is approximately 1 Hz at an epoch 
length of approximately 1 second, there are more than 150 spectral lines 


lying between two harmonics whereby said spectral lines contain only the 
noise resulting from the asynchronous EEC If one applies the statistical 
test to these spectral lines, then one can obtain approximately 3,000,000 
test values at 100 recordings of each 200 epochs (100x200x150). The 
distribution computed from these test values represent a very good 
estimate of the density function of the null hypothesis fitting the real data 
from which the searched critical test value can be read. 

The density function of the null hypothesis calculated by means of 
the hereby described method is shown in Figure 2 as an example for the 
Test Modification 1. 

The use of the inventive solution has the following advantage: 

While the known solution uses only partially the information 
contained in the spectrum because of the limitation of the ranks of the 
phase angles (q-sample uniform scores test) or the ranks of the phase 
angles together with the spectral amplitudes (Modification 1), more 
information (Modification 2 and 3) or the entire information contained in 
the spectrum- (Modification 4) is used with the inventive modifications 
(Modification 2-4). 

The result is a higher test power of the proposed modifications. A 
higher test power leads to the fact that the responses are detected more 
rapidly during hearing screening with a preset stimulus level and the time 
for screening is thereby shorter. A more exact objective threshold 
determination is made possible in an objective audible threshold 
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determination as a result of the higher test power since the response 
verification becomes closer to the audible threshold of the patient himself. 

The elimination of the distribution independence in the proposed 
modification is not a disadvantage since the appropriate distribution of the 
null hypothesis is determined by the data specified in the proposed 
method. 


